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ETT DATA PROCESSING SYl 

This invention relates to the field of electrical impedance tomography (ETT), and in particular to 
ETT data processing systems and methods. 

There are many industrial two-component flow applications where it is required to monitor the 
flow rate of one, or both, of the flowing components, la a two component flow where the 
continuous component is electrically conductive, and where there is a conductivity contrast 
between the continuous and discontinuous components, a technique known as dual-plane Electrical 
Impedance Tomography (ETT) has the potential to be used on-line to obtain accurate estimates of 
the local disperse phase volumetric flow rate, the mean disperse and continuous phase volume 
tractions and the distributions of the local axial, radial and angular velocity components of the 



IS The invention is concerned with online data acquisition and processing of measurements from two- 
phase flows* whose axial velocities can be up to 10ms" 1 , with ±5% precision. To be able to measure 
two-phase flows at these velocities requires data acquisition and processing speeds of 1000 frames 
per second per dual-plane in order to generate flow velocity maps and this has not been possible to 
date, using conventional techniques. 

20 * 

The following prior art documents are incorporated herein by reference. 

DICKIN F, WANG M, Electrical resistance tomography for process applications, Meas. ScL 
Technol. 7 (1990, 247-260 

WANG M, DICKIN FJ, BECK MS, Improved electrical impedance tomography data collection 
system and -measurement protocols, Tomography Technique and Process Design and Operation 
1SBN1813 122467 

30 TTI, (2000), TMS320C6202/B Fixed-Poira Digital Signal Processors, Texas Instruments 
Incorporated 



WANG M, (1994), Electrical impedance tomography on conducting walled process vessels, PhD 
thesis, UMIST 

SUMMARY OF THE INVENTION 

It is an object of the present invention to provide systems and methods for multiphase flow 
measuring that are robust, non-invasive and can deliver on-line multi-dimensional flow 
information- 



According to a first aspect of the invention, there is provided an electrical impedance tomography 
(ETT) data processing system, for acquiring and processing data from two-phase flows, comprising 
a dual-plane sensor, a plurality of digital signal processing modules configured in a data pipeline 
45 processing ar r an gement and a plurality of data acquisition subsystems in communication with a 
first one of said digital signal processing modules. 
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Preferably;, the data pipeline processing arrangement is capable of acquiring and/or processing 
1000 < 



5 Preferably, said dual-plan© sensor comprises an electrode array in communication with one of said 
data acquisition subsystems. 

In a preferred embodiment, said data acquisition subsystems each include one or more of a voltage 
controlled current source, an equal-width pulse synthesiser, a synchronised digital demodulation 
10 unit and an over-zero switch. 

The voltage controlled current source preferably comprises a parallel structure of eight AD 8 44 
chips or equivalents, configured as four pairs in which the two inverting inputs of each pair are 
cascaded together with a current-setting resistor. The four negative and four positive current. 
15 outputs of the AD844 chips are summed together to provide total negative and positive c u ii ent 
outputs. 

Preferably, the voltage controlled current source further comprises a DC restore fecility in which a 
capacitor and a resistor connected to the non-inverting input of each AD 844 chip restore DC 
20 components and cancel Hie DC offset at the current outputs. 

In a preferred embodiment, the voltage controlled current source includes a potentiometer for 
- amplitude balancing. 

25 Preferably, said equal-width pulse synthesiser comprises a clock signal for triggering an address 
generator to continuously output addresses to a pre-programmed memory, the memory output 
being connected to a digital to analogue converter, the digital to analogue converter providing a 
staircase signal output, characterised in that different sampling rates are provided at different sign al 
frequencies so that the staircase signal output has the same time step-length at all frequencies. 

30 ' 

Preferably, the synchronised digital demodulation unit comprises 16 sets of programmable gain 
amplifiers each having an analogue to digital converter, a strobed First-In First-Out (FIFO) 
memory and control logic and "wherein the signal output of the equal-width pulse synthesiser 
triggers said 16 analogue to digital converters to acquire measurement data in parallel at predefined 

35 intervals. 

Preferably, said over-zero switch comprises two multiplexers, each having 16 electrodes mounted 
in one sensing plane connected to their 16 output channels and two flip-flops controlled by a 
switching command from one of said DSP modules and, ideally, a D-type flip-flap for preventing 
40 data being written prematurely to the second of said flip-flops. 
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According to a third aspect of the invention there i$ provided a recording medium having recorded 
thereon, computer implementable instructions for performing the method described in the preceding, 
paragraph* 

The system has been carefully designed in a modular fashion and can consist of several data 
acquisition modules and computing modules. The data acquisition modules include a voltage 
controlled current source with a DC-restore circuit, an equal-width pulse synthesizer unit to 
produce synthetic waveforms for electric field excitation, synchronized switchings and sampling 
control, and a synchronized digital demodulation unit consisting of sixteen parallel measurement 
channels. The computing modules include a powerful digital signal processor (for example 
XMS320C6202B/671S) with an 3EEE1394 communication interface. Several DSP modules can be 
pipelined for a series of tasks ranging from measurement control to image reconstruction to flow 
velocity implementation. The DSP system can implement two image reconstruction methods: the 
single-step SCG method and the back-projection method. An online-updating crossrcorrelafion 
algorithm has lateen developed and is to be implemented on the system to calculate flow velocity. 
This technique is less time and memory consuming compared to direct cross-correlation 
techniques. The system architecture, detailed modular design, and data processing algorithms are 
, introduced below. With a litrle optimization on hardware, a speed of data acquisition of 1 1 64 dual- 
frames (2383 million data points) per second has been achieved with a root mean square value less 
^ than 0.8%@80kHz using a water phantom. 

. The very high speed of data acquisition and processing afforded by the system of the present 
invention is achieved as a result of one or more of the following: 
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reducing the time required for multiple-tasks, such as data collection, image reconstruction 
and flow velocity implementation by the use of the data pipeline process with multi DSP 
(see Figure 1); 



• the speed of signal coupling due to the transient time of the improved wide bandwidth 
coupling circuit with differential discharge circuits, therefore significantly reducing the 
transient time (see the SDD, Figure 10); 

• the residual voltage in capacitance of coupling due to the electric charge being switched at 
~ 'the zero voltage point; or zero charged energy of the capacitance charged by the sinusoidal 
^ current, significantly reducing the transient time (shown in the OZS of Figure 1 1); 



the data collection operation speed, improved by the parallel structure of the SDD of 
Figure 8; 

the time required for transferring and storing data, to memory, reduced by the logical 
controlled FEFO memory (see the SDD and Figures 3 and 8); 



• the limited sampling speed of the analogue to digital converter T hrift? the application to 
^ high frequency in conventional methods but this limitation was overcome in the present 

system by employing a rolling sampling method (see SDD and Figure 9d); 
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• the time limitation of the control operation of the micro-processor, improved by the use of 
the logical control (not DSP), see SDD and Figure 8; 

• improved data transfer from data acquisition subsystems to DSP to empty FIFO memory 
5 for the next operation, realized by the use of a fast DSP-FIFO interfere (Figure I). 

Using the system and method of the present invention, properties of the flew can be measured 
accurately even when the flow is highly noiHinifbrm. Dual-plane EXT has the added advantage of 
providing on-line volume fraction and velocity distributions which can be used for the purposes of 
10 process monitoring and control. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Preferred embodiments of the present invention will now be more particularly described,, by way of 
15 example only, with reference to the accompanying drawings, in which: 

Figure 1 illustrates the system architecture; 

Figure 2 shows one possible system configuration of a high performance two-phase flow meter 
20 using the DSP module as a building block, (a) Single Read-Operation Timing, (b) Single Write- 
Operation Timing; 

Figure 3 illustrates the architecture of the BIT data acquisition subsystems (DAS); 
25 Figure 4 shows a simplified schematic (a) and an equivalent schematic (b) for AD844; 
Figure 5 shows a voltage controlled current source (VCCS); 
Figure 6 shows the structure of the Equal-Width Pulse Synthesizer (EWPS); 

30 

Figure 7 shows waveforms generated by the EWPS; 

Figure 8 shows the structure of the Synchronised Digital Demodulation (SDD); 

35 Figure 9 shows the configuration of the sampling modes, (a) 4 samplings, (b) 8 sam plings , (c) 16 
samplings, (d) rolling-sample method; 

Figure 10 shows an AC-coupling network; 

40 Figure 11 shows the trigger circuit of the over zero switch (OZS); 

Figure 12 shows a typical curve of independent measurements obtained from a phantom filled with 
mains tap water; 

45 Figure 13(a) shows photographs of two balls dropping, erne relatively heavy and one relatively 
light; * ~ 
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Figure 13(b) shows 2D ECT images; 

Figure 13(c). show 1000 stacked images along section A in Fig 13(b) where a relatively light ball 
5 is falling through; 

Figure 13(d). shows 1000 stacked images along section B in Fig 13(b) where a relative^ heavy ball 
is felling through; 

10 Figure 13(e) shows the variation of averaged conductivities at Hie pixels along the sections A and 
B, according to time elapsed; 
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Figure 14 shpws cross-correlation curves at fits; sampling pixels along Hie radius of sensing planes 
via data collection speeds; 

Figure 15 shows 2D flow velocity profiles from an oil-water axial flow (M4) with superficial 
velocity: Vwater=0^49m/s, Voil= 0.05m/s; and 



Figure 16 shows flow velocity profiles from an air-water swirling flow with superficial velocity; 
20 Vwater=0.39m/s, Vgas=0»067m/s. (a) 2D profile from the combination of angular and radial 
velocity; (b) 3D profiles from the combination of all 3 dimensional velocity profiles. 



DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Throughout these documents, the following terms are used: 
"on-line measurement" — synchronised measuring and imaging on the process line* 

30 

"two-phase flow" - a fluid comprising two components, wherein there are differences in the 
electrical properties of the two components. 

W SCG method 33 — a sensitivity theorem based algorithm using a conjugate grading method. 

3 S "back projection method 33 — a linear approximation algorithm. 

"data pipeline processing arrangement 33 — an arrangement fax which each processor (DSP) in the 
pipeline performs a particular function on the data and then hands over to the next DSP in the 
pipeline. All of the processors operate on data simultaneously so that data can be processed much 
40 more quickly than would be possible with a conventional DSP arrangement. 

System Overview 

The measurement system, as shown in Figure 1, includes up to four Digital Signal Processor 
45 modules (DSP1-DSP4), two identical Data Acquisition Subsystems (DAS1, DAS2) connected to 
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two electrode arrays (Electrode Array 1, Electrode Array 2), and a third Data. Acquisition 
- Subsystem (DAS3) collecting data from conductivity, pressure, temperature and other auxiliary 
sensors/transducers. All three DASs are controlled by system-controlling DSP module DSP1 
through a common interface (Data Acquisition Interface). Each of DSP1-DSP4 can be configured 
5 to communicate one at a time with a remote PC, 

All of the DSP modules are preferably identical to reduce the design and debugging tasks. The 
1MS320C6202B/6713 DSP chip [IT 2000] may be used for the DSP module. Any DSP module 
can communicate with a remote PC via an IEEE 1394 interface, shown in Figure 1, to make the 
10 functionality of the system flexible. The four DSP modules work on a data pipeline processing 
technique to increase the data^processing speed, 

They can perform different tasks on data in the pipeline, for example in the illustrated embodiment, 
DSP1 controls data acquisition and processes raw data; DSP2 performs image reconstruction, 
15 DSP3 fuses the image data to obtain flow information by performing velocity calculations and 
DSP4 is provided to carry out other specified functions. 

Figure 2 is a block diagram showing how DSP modules are used to implement such a high 
performance system. DSP1 acts as the system controller and controls the data acquisition process 
20 of the three data acquisition subsystems (DAS1-DAS3) which are only connected to DSPL DSP1 
then formats the data and transfers it to DSP2. The data are transmitted via McBSPO on DSP 1 and 
received via McBSPl on DSP2. Besides these tasks, DSP1 is also responsible for the IEBE1394 
interface communication! 

25 EFT data acquisition subsystems 

Two identical d#a acquisition subsystems (DAS1, DAS2) are included in the system to collect data 
from two electrode arrays (a dual-plane sensor). A block diagram of a DAS module is shown inside 
the dashed-line box in Figure 3* The four modules of most importance in the DAS module are: a 
30 Voltage Controlled Current Source (VCCS), an Equal-Width Pulse Synthesizer. (EWPS), a 
Synchronized Digital Demodulation (SDD) unit and an over-zero switch (OZS) in the control base. 

Voltage Controlled Current Source (VCCS) 

35 The AD 844 from Analogue Devices has a wide bandwidth of 60 MHz and very high output slew 
rate of up to 2000 V/ms. It differs from a conventional op amp in that the signal inputs have 
radically different impedances. As schematics shown in Figure 4, its non-inverting input (Pin 3) 
presents the unusual high impedance while the inverting input (Pin 2) presents low impedance, R©j, 
of about S0SX A current applied to the inverting input is transferred to a complementary pair of 

4Q wideband unity-gain current mirrors which deliver the same current to an internal slewing (TZ) 
node (Pin 5). It can be seen that Pin 5 and the inputs behave as a good VCCS with high output 
impedance and wideband. The one drawback is its small transient current range of ±10mA in 
comparison with general ERT requirements. To overcome this problem, a parallel structure 
consisting of eight AD844s has been developed and shown in Figure 5, 

45 
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Two inverse-phased sin-wave voltage inputs are buffered and then connected to the non-inverting 
input of each AD844, where two capacitances C and two resistances r can perfectly restore the 
Direct-Current (DC) components generated in the circuit and then cancel the DC offset at the 
current outputs* Eight AD844s constitute 4 pairs of two-output VCCSs. Two inverting inputs of 
5 each pair of AD 844s are cascaded together with a enrcent-setting resistor R. The positive and 
negative current outputs of each pair are summed together to form larger current outputs. The total 
output currents can be estimated using the following function: 

±IauT s =±2^V SIN /(2xR 1Kf +R) (Equation 1) 

10 

Compared to other VCCS designs, the advantage of the VCCS shown in Figure 5 is its excellent 
balance property. The ±1^ have ideal 180° phase difference* and ignorable DC biasing and then* 

amplitudes can be easily balanced with the potentiometer;, Rp. The maximum output current 
amplitude of this circuit is 30mA. Its output impedance is about 750k £2 []18pF. 

15 

Equal-Width Pulse Synthesizer (EWPS) 

The EWPS structure, shown in Figure 6, is similar to that of the widely used Equal-Sample Signal 
Generator (ESSG). An auto-run clock signal triggers the address generator to continuously output 
20 addresses to a pre-programmed memoiy. The memory then sequentially outputs data to a 12-bit 
Digital-to-Analogue Converter (DAC). A low-pass filter then smoothes the staircase signal output 
from the DAC to the required sine-wave voltage signal. The frequency, phase and amplitude of the 
generated signal can be set by outside controlling unit DSP. 

25 The difference between EWPS and ESSG is that the former has different saniple numbers at 
different frequencies to make the staircase wave use the same time step-length at all frequencies, 
while the latter has the same number of samples at all frequencies so the staircase wave has 
different time step-lengths at different frequencies. This difference has led to a significant reduction 
in the complexity of the filter design using EWPS as it only needs to set one comer-frequency for 

30 the filter for all signal frequencies while the ESSG needs to set different corner-frequencies for the 
filter for different signal frequencies. 

Considering the use of the same step-length, the sampling pake frequency f 9 is unchanged 

/, or / 0 =Lf s (Equation 2) 

35 Therefore, 

^^^M^^Wf) ' (Equations) 

Hence, the lowest harmonics from the EWPS is N ' and the EWPS may use one low pass filter 
with a fixed comer frequency for all waveforms. The step-length used in the EWPS is 0*5 fis which 
40 gives a sampling pulse frequency of 2.0MHz. The waveforms generated by EWPS are illustrated in 
Figure 7. 



7 



WO 2005/022138 



PCT/GB2004/003670 



The specialty designed EWPS also provides a convenient way to produce synchronized signals and 
sample control signals that are necessary in controlling SDD. Hie synchronous signal (square wave) 
is used to control die start and the end of the sampling sequences and the sampling pulses are used 
to actually trigger ADC*s conversions at specific times. The timing of both signals and the width of 
5 sampling pulses play a very important role in the achievement of SDD accuracy. 

Synchronized Digital Demodulator (SDD) 

The Synchronized Digital Demodulator (SDD) is a high speed, synchronized parallel digital 
10 demodulation system, which is composed of 16 sets of Programmable Gain Amplifiers (PGA) and 
12-bit ADC, two sets of Strobed First-Bi First-Out (FIFO) Memory and control logic. Figure 8 
shows the functional block diagram of the SDD. The sampling pulse signal fiom EWPS triggers all 
sixteen ADCs to acquire measurement data in parallel at precisely defined intervals. Two sampling 
modes are available. 

15 

In Mode I, 4, 8, 16, or 32 samples can be selected in one sine, wave 4 cycle as the same as the 
conventional sampling method. ' 

Mode II uses a rolling-sample method. It acquires all samples over several sine-wave periods with 
20 one sample per period and a fixed phase shift between one sample pulse and the one after. 
Sampling Mode H is particularly designed for sampling; at a high frequency due to the speed 
limitation of the ADC. 

The configuration of the two sampling modes is shown in Figure 9. The completion signal 
25 combined from ADC's conversion status signals will then trigger a control clock to sequentially 
push converted digital data into two FIFO memories. All operations are fully logical controlled 
without the involvement of the DSP, except the starting signal of each frame of measurement 

The coupling circuits, as shown in Figure 10, are made of a number of voltage buffers and passive 
30 components to provide a wide bandwidth and fast setding ac-couplmg interface, as well as a do 
bypass to the dc current bias of the programmable gain amplifier (PGA). lie input can' be 
configured in either the single or the differential input mode. Only the differential mode is 
illustrated in Figure 10 for the use of the adjacent measurement strategy in electrical impedance 
tomography, where R1-R3 and Rl is much smaller than K2. Therefore, R2, CI dominate the 
35 bandwidth of the coupling, Rl and R3 provide a DC bypass to the DC current offset of PGA but 
the responding speed is fast since the transient time is dominated by Rl oxR3. 

Over Zero Switch (OZS) 

40 Conventional ERT systems use sinusoidal current sources to excite the sensing field. Since a 
controlling unit, such as a PC, DSP or microcomputer, directly controls the current switching 
between electrodes, the switching operation may occur randomly at any time during a sinusoidal 
cycle. When the current value at the time of the switching operation is non-zero, electrical charge 
may reside in the capacitances of the signal passing circuit, such as the double layer capacitance at 

45 the electrode-electrolyte interface, the coupling capacitance and the stray capacitance, where the 
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double capacitance plays a major role in all capacitances. These may result in residual potentials 
presented between these capacitances and produce an additional contribution to next excitation and 
measurement These residual charges will require several sinusoidal cycles to die-out, which means 
the whole data acquisition system (PAS) has to wait for several cycles after each switching 
S operation in order to perform next voltage measurement, 

The over-zero switch (OZS) aims to eliminate the charged residual voltage by controlling the 
switching operation of the current inj ection. In an ideal cases no charged potential resides in either 
the double layer capacitance or the coupling capacitance after switching off the current injection so 

10 ^ settling time is shorter with the data acquisition speed being increased. Since the charged 
voltage has a it/2 delay with respect to the current, the switch should operated at the time either a 
zero-voltage presents between the double-layer capacitanoe or a nmimum-crareiit through the 
capacitance. The ourrent-marimnxn is actually applied to optimise the switching time or the 
relevant phase difference between the synchronised signal and the sinusoidal signal, as the charged 

15 voltage between fee double-layer capacitance is hardly measurable. 

Figure 11 is the circuit used in the improved system to realize the over- current-zero switching 
technique. Two inverse-phased currents 1+ & I- are connected to the common inputs of two 
multiplexers Ul & U2 respectively, 16 electrodes mounted in one sensing plane are connected to 

20 the 16 output channels (CH[0:15]) of each multiplexer* Channel selection is together controlled by 
DSP command and a synchronize signal, SYN_SIK. SYNJS1N is a square-wave signal, which has 
the same over-zero point with two current signals 1+ & I-- When DSP send out a switching 
command, an 8-bits data (each multiplexer is controlled by 4-bits data), ED[0;7L is written into the 
first flip-flop, U3. This data is later written into the second flip-flop, U4, by the next rising edge of 

25 SYNJSN, and then make fee two multiplexers switch to other channels. A D-type flip-flop, U5A, 
is used in this circuit to prevent writing data into U4 before the data has been set-up at the output of 
IB. 



30 . Test results 

A laboratory test vessel (a "phantom 7 ') fitted with 16 stainless steel electrodes and filled with mains 
tap water was used in one of the trials. In the experiment, all the PGAs in the SDD were equally set 
to a fixed gain value in order to simplify the test The current excitation frequency was 80 kHz. 
3 5 Measurements were taken from all electrodes during each current excitation/projection, which gave 
256 measurements. Four samples per period were taken. Therefore, each frame data has a total of 
1024 measurements. A simplified digital matched-filter was applied to demodulate the data. A set 
of independent measurements were then extracted and averaged from their mirror data based on the 
reciprocal theorem. 

40 

The top and middle curves in Kgwe 12 present two of measurements obtained from vessel filled 
with mains tap water. The bottom curve denotes the relative change between the two measurements 
to reflect the repeatability of these measurements, To demonstrate the high speed of the new ETT 
system, a test was earned out by dropping two non-conductive balls. The two balls have similar 
45 dimensions but difierent densities (radius, SG and conductivity of the two balls are as follows: 
Rl^L74cm, R2=L75cm, S(H=2*55 9 SG2=1.13 and o\ 2=0). The two balls were dropped to the 
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bottom of the phantom (see Fig. 13a), arriving at different times- The time difference can be 
• dearly visualized from the HT images (Fig>13b), 1000 stacked ETT images along the sections A 
and B (Fig.l3c and 13d) and the reconstructed conductivity changes at relevant pixels as shown in 
Fig.l3e. The relatively heavy ball dropped in section A arrived at the bottom earlier than the 
5 relatively light ball dropped at section B, It can even be seen from both the photographs and ETT 
images (Figs.l3d & e) that Ac heavy ball bounced back when it first knocked the bottom. 

The total n umb er of data collection points depends on the sampling rate (the number of samples per 
signal period). With a fixed sampling rate of 320kHz, there are 1024 measurement points from on© 

10 sensing plane with 16 electrodes for one frame image at the signal frequency 80kHz, or 2.048 
million data collection points from a dual-plane sensor for 2000 frames of images at the signal 
frequency 80kHz, After the optimization of the controlling unit of DSP C6202 modules at aspects 
of data collection procedure, memory management; internal data transferring and program coding, 
die data acquisition speed of the system increases significantly. A data acquisition speed of 1163.6 

IS dual-frames per second (dips) at 80kHz (23212 frames per second in total) can be achieved. The 
data acquisition speeds at d ifferen t signal frequencies are summarized below: 
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Data Processing Algorithms 

Two image reconstruction algorithms have been coded for the DSP module. They are Sensitrvity- 
map-based Bade Projection (SBP) method and a simplified SCG method - Single-step SCO 
(SSCG) method. The SSCG method is based on the SCG method but only iterated once from a 
homogenous conductivity distribution. The SBP method is the simplest and fastest It seems that 
the SSCG has a better performance than the SBP method. 

Data format in floating point is adopted in most cases due to the use of the ill-conditioned 
sensitivity matrix for image reconstruction. However, a DSP module C6202 with fixed-point data 
format was used in an earlier version ofthe system which was relatively slow to process the data in 
floating-point format. It only can reconstruct images with the SBP method at a speed of 189.0 
frames per second, which is much lower than the online data collection speed. Therefore* a 
floating-point DSP module, C6713 is used in the system ofthe present invention to obtain a higher 
image reconstruction speed. This has achieved an image reconstruction rate of 577.4 frames per 
second 



An onlme^updating cross-correlation algorithm has teen developed Assuming that the two ETT 
electrode arrays are installed parallel on the flow-pipe, the axial-flow velocity distribution can be 
estimated by a direct cross-correlation method, as given in Equation 4: 

Az(»)- «-<W.~.<N-» (Equation 4) 
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15 



where N is the sample length, and fi(m) & fe(m) are the m* up-flow and down-flow images 
respectively. Eq.l can be simply implemented online by updating the Riz(n) with the new k-th 
images, as described in the following: 

where superscript (k) indicates thatfbe corresponding value is of k-th sample, This implementation 
can greatly save calculation time and reduce memory size. The computation speed is subject to the 
number of pairs of images that are estimated as a maximum sampling length. 

The minimum acquisition time, t, can be taken as twice Hie product of the minimum transit time of 
the fluid, and the fractional velocity discrimination, k. Rearrangement gives the following 

k-t/2 (Equation 6) 



Therefore, the speed of data collection plays a dominant role in the use of cross-couelation 
algorithm. If the distance between two sensing planes is d — 0.1m then O.lm/lOms^ 1 — 0.01s. Iff 
= 1000 frames/s was achieved then t - 0.001s. Using Equation 6 gives a discriminatory precision of 
0.05 or 5%. The local velocity can be derived from the correlated frame number, n, at each im agi n g 
^0 pixel, as given by Equation 7: 

c = fd/n (Equation 7) 



The importance of data collection speed is demonstrated in Figure 14. The legend numbers are the 
x- and y- coordinates for a vessel with a diameter 20. Due to the limited data collection speed 
(Figure 14(a)), the correlated frame numbers (in respect to the lowest cc-value) are small and 
therefore, they produce a poor resolution in flow velocity measurement. A good resolution from the 
fast system is clearly illustrated in Figure 14(b). 
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